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Abstract A quasi-non-linear fracture mechanics model 
based on beam on elastic foundation theory is applied for 
analysis of the splitting failure of dowel joints loaded per¬ 
pendicular to grain. Simply supported beams symmetrically 
loaded by two dowels are considered, and the effects of 
edge distance, dowel spacing, and distance between dowels 
and supports are accounted for. The foundation modulus 
used in the beam on elastic foundation model is chosen so 
that the perpendicular-to-grain tensile strength and fracture 
energy properties of the wood are correctly represented. 
This ensures that a conventional stress analysis and failure 
criterion lead to the same solution as the compliance 
method of fracture mechanics. A semiempirical efficiency 
factor is proposed to account for the influence of the total 
beam depth, which does not enter the beam on elastic foun¬ 
dation model, but the effect of which is evident from tests. 
It is shown that the so-called Van der Put/Leijten model, 
which recently has been adopted in Eurocode 5, appears 
as a special case of the model presented. Tests on simply 
supported beams with a single dowel joint at midspan are 
compared with the theoretical predictions. Various edge 
distances, beam depths, and spans were tested. 

Key words Dowel joints ■ Perpendicular to grain • Split¬ 
ting ■ Fracture mechanics ■ Beam on elastic foundation 


Introduction 

A dowel-type fastener joint loading a beam perpendicular 
to grain may fail in a ductile manner, characterized by bend- 
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ing of the fasteners and/or embedment of the fasteners into 
the wood, or it may cause a brittle failure in the beam 
characterized by splitting of the wood. The ductile failure 
modes are well understood and can be predicted fairly accu¬ 
rately by the European Yield Model, 1 which now forms the 
basis of design of dowel-type fastener joints in major design 
codes. Only recently, a number of simple analytical models 
based on fracture mechanics have been proposed 2-7 for 
prediction of brittle failure, but none has yet gained wide 
acceptance. 

A quasi-non-linear fracture mechanics model based on 
beam on elastic foundation (BEF) theory was previously 
presented. 6 The model presented in this article is based on 
the same ideas, but is extended to take into account the 
influence of the supports in simply supported beams. 

Application of the BEF theory can only be expected to 
lead to reasonable results for edge distances, h e , which are 
relatively small compared with the total beam depth, h. 
By nature, the BEF model does not predict any influence of 
h. However, tests clearly indicate the counterintuitive effect 
that the failure load increases if h decreases for constant h c . 
In the present article a semiempirical approach is suggested 
for modifying the model to produce reasonable solutions 
for larger hjh ratios (a). The suggested modification makes 
the presented model contain previous models 2,5 as special 
cases. 


Theory 

A structure as shown in Fig. 1 is considered. The part of the 
beam below the dowels in Fig. 1 is supposed to behave as a 
beam (beam depth: A e ) on an elastic foundation. The part 
of the beam above the dowels is supposed to provide the 
foundation. 

The symmetrical half of the structure is modeled as 
shown in Fig. 2. Note that the structure analyzed previously 6 
(hereafter denoted as “plate joint"’) differs from Fig. 2 by 
having no point support at the right end. 
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Fig. 1. Geometry of structure 



The foundation stress, o(s ) (tensile stresses are positive), 
at the location of the dowel is given by 


Fig. 2. Beam on elastic foundation model 


cr(s) = — Kw(s) 


(1) Case 2: 


where K (unit: N/m 3 ) is the foundation modulus and vv(.v) is 
the deflection of the beam axis. 8 

Failure is assumed to occur when the stress, o(s), equals 
the tensile strength perpendicular to grain, /„ of the wood, 
leading to the failure load, P c . 

The foundation is intended to model the strength and 
fracture performance of the wood perpendicular to grain. 
The damage and fracture performance of wood is non¬ 
linear, but is in the present analysis represented by a linear 
response that is equivalent in terms of peak stress, f , and 
fracture energy dissipation, (g t . Because the energy dissipa¬ 
tion in the case of linear performance is 

il follows that K = ( 2 ) 

The solution to the governing differential equations for a 
Timoshenko beam on elastic foundation falls into two cases, 
the limit between which may be characterized by a require¬ 
ment to the depth, h„ of the beam on elastic foundation 
(rectangular cross section) 


Case 1: h s 


200 G 2 g t 
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Case 2: h. 
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where E is the modulus of elasticity in the grain direction 
and G is the shear modulus. 

For the beam shown in Fig. 2, the failure load is found to 
be 


Case 1: 
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b being the width of the beam, I the moment of inertia, and 
A s is the equivalent shear area of the beam cross section 
(A, = 5bhJ 6 and I = bh 3 J 12 for a rectangular cross section 
with beam depth hf). 


Infinitely long beam with two dowels 


For an infinitely long beam (L —» °°), the failure load is 
reduced to (for u u u 2 , v u and v 2 , see Eq. 4) 
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Case 2: P c = 2 bf - 


vl(l + e^ s ) - u 3 2 ( 1 + e -2 " 2 *) 
which is the same solution as previously obtained. 6 
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Fig. 3. Plate joint and beam joint solutions 


Fig. 4. Effect of dowel spacing and supports 


Finite beam with one dowel 


The failure load, P c , for a beam with a single dowel is 
obtained from Eq. 4 by setting s = 0 (and multiplying by 2). 
For u and v, see Eq. 4. 


Case 1: P c = 8 bf t 


Case 2: P c = 2 bf t 
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Infinitely long beam with one dowel 


Figure 3 shows the investigation of a specimen with a 
single dowel (or joint) at midspan. The beam joint solution 
(Eq. 7) is compared with the plate joint solution 6 (i.e., no 
point support at the right end in Fig. 2). Examples showing 
case 1 and case 2 solutions are given. In both examples, E = 
12700MPa, G = 870Mpa, and edge distance h c = 56mm. 
For case 1, the fracture properties reported for laminated 
veneer lumber (LVL), 5 / t = 0.89MPa, (g t = 0.25 N/mm, are 
used. For case 2, f = 5.0MPa and Q t = 0.30N/mm. 

In Fig. 4, a simply supported beam with two dowels (or 
joints) is examined. The effect of dowel spacing, 2s, and the 
influence of supports according to Eq. 4 is shown for a beam 
of length L = 15h e (h,. = 56 mm) using the same material 
properties as in Fig. 3. 


For an infinitely long beam with a single dowel, the previ¬ 
ously obtained solution 6 is arrived at 


Pq — yPc,LE ™ 


( 8 ) 


where the linear elastic fracture mechanics solution, / j c L ef M , 
is given by 


-fc.LEFM — 2 bC 12 Jh^, 


and 7 is given by 
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( 10 ) 


The above failure loads were derived using a conven¬ 
tional stress analysis and failure criterion, and fracture me¬ 
chanics enters the model through the constitutive relation 
adopted for the foundation. This analysis is an exact anal¬ 
ogy to the quasi-non-linear analysis of lap joints . 9 

It has previously been demonstrated that the linearized 
constitutive assumption (Eq. 2) used in the BEF model 
makes the stress approach lead to the same solution 
obtained by the compliance method . 10,11 


Modification of foundation stiffness 


In the BEF theory described above, the foundation is as¬ 
signed the stiffness as given by Eq. 2. This may be a good 
approximation for small values of a = hjh , but for large 
values of a the completely stiff support of the springs as 
shown in Fig. 2 is not a realistic assumption. 

K s denotes an additional contribution to the foundation 
stiffness due to the deflection of the part of the beam above 
the potential crack path (beam depth: (1 — a)h). It is here 
assumed that K s can be written as 


K, = K\ 


1 — a 


( 11 ) 


where a = hjh , K is given by Eq. 2, and m is a parameter to 
be determined empirically. Note that m = 1 corresponds to 
KJK = AJA s2 , and m = 3 corresponds to KJK = EIJEI 2 , 
where subscript 1 refers to the beam part above the poten¬ 
tial crack path (foundation) and subscript 2 refers to the 
beam part below the potential crack path (beam on elastic 
foundation). 

The modified foundation stiffness is given by 
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Fig. 5. Effect of modification of foundation stiffness 
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Setting m = 1 results in (p = 1 — a, i.e., 
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The additional foundation stiffness as introduced by Eq. 11 
is constant along the beam axis, whereas a proper account 
for the additional deflections of the part of the beam above 
the crack (beam depth: (1 — a)h) would involve a founda¬ 
tion stiffness that varies along the beam axis. The 
modification of the foundation stiffness as suggested here 
may therefore be considered a semiempirical attempt to 
involve the influence of the total beam depth, h. 

Figure 5 shows an example of the influence of the modi¬ 
fied foundation stiffness as given by Eq. 14. A simply 
supported beam with a single joint at midspan (b = 25 mm, 
h s = 40 mm, E = 7880 MPa, G = 438 MPa,/, = 1.2 MPa, and 
(g t = 0.16N/mm) was analyzed by means of Eq. 7. The 
efficiency factor, k, is defined as k = P c (K moi )/PJK), where 
PJK mo J) is the failure load obtained using K mod and PJK) is 
the failure load using K. 

Figure 5 shows that the efficiency factor depends not 
only on the a value ( hjh ), but also on the L/h c ratio. For 
Llh c —> °°, k asymptotically approaches k, . 
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Instead of using the modified value of the fracture en¬ 
ergy as given by Eq. 14, the modification factor, tc„, as given 
by Eq. 16, may be used as a semiempirical efficiency factor, 
by which the solutions given by Eqs. 4, 6 , 7, and 8 may be 
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Fig. 6. Test setup 

multiplied to account for the effect of the total beam depth, 
h. Although theoretically not strictly conservative for all U 
h e values (e.g., according to Fig. 5, k assumes a minimum 
value for L//z e values between 4 and 6 for a = 0.8), k„ seems 
to be a practical and safe simplification. 

It is noted that —> p for f —> 00 . The Van der 

Put/Leijten solution 2,7 is thus obtained by multiplying Eq. 8 
by K,_ and assuming/, —» 00 . The Van der Put/Leijten solution 
seems to give a prediction of the influence of the total 
beam depth, h, in good agreement with test results . 2 ^ 1 It 
is further noted that —> 1 for a = hjh —» 0. The 

Gustafsson/Larsen model 5 is thus obtained by assuming 
a —> 0 and f —> °°. 


Experimental 

Materials and methods 

Simply supported glulam beams of Japanese cedar ( Cryp- 
tomeria japonica ) were tested using joints with a single 
dowel as shown in Fig. 6 . 

The glulam was made of 30-mm laminae of the same 
grade (machine graded). Specimens were cut from eight 
beams (100 X 200mm 2 cross section) of which the modulus 
of elasticity (MOE) was measured (longitudinal vibration 
frequency). MOE ranged from 7500MPa to 8250MPa 
(mean value: 7880MPa). Moisture content (MC) was 
12.5%. Density ranged from 362kg/m" to 392 kg/m" (mean 
value: 373 kg/m 3 ) at the given MC. Only knot-free speci¬ 
mens were used in order to limit the variation. In all cases 
beam width was b = 25 mm, and a 14-mm dowel was used in 
a 15-mm hole. The load was displacement controlled, and 
time to failure was 2-3 min. 

Two different edge distances, h,. = 20 mm and h e = 
40 mm, were tested. Edge distance h c = 20 mm was tested 
for h = 200 mm and h — 100 mm (i.e., a = hjh = 0.1 and 
0.2). For each beam depth L//z e = 2.5, 5, 7.5, and 10 were 
tested. Edge distance h : . = 40mm was tested for h = 200mm 
and h = 80mm (i.e., a = hjh — 0.2 and 0.5). For 
h = 200mm, Llh e = 2.5, 5, and 23.75 were tested. For 
h = 80mm, Llh e = 2.5, 5, and 7.5 were tested. Nine to 12 
specimens were tested for each condition. 

Tests were also conducted on double cantilever beam 
(DCB) specimens and on so-called plate-joint specimens, as 
shown in Fig. 7. All specimens were selected to avoid knots. 
The DCB specimens were used to determine the fracture 
energy. Thirteen specimens were tested, b = 25 mm, 
h = 40 mm, L = 500 mm, a = 200 mm. 

The plate joint specimens were tested for comparison 
with the beam tests ( b = 25mm, h = 200mm, L = 550mm). 
Eight specimens were tested for h e = 20 mm and /z c = 40 mm. 
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Fig. 7. Geometry of double cantilever beam (a) and plate joint 
(b) specimens 


Results and discussion 

Failure load 

Fracture energy was determined from the DCB test results 
using Eq. 17,' 5 resulting in (g t = 0.16 ± 0.05 N/mm. 

G, = 12-A—(« + daf, da = hj— (17) 

y£ b 2 h 3 E [ ’ AllO G ' 

The plate joint tests resulted in: P c = 2054 ± 336N for h c = 
20mm, and P c = 2867 ± 371N for /t e = 40mm. 

The results of the beam tests are shown in Fig. 8 along 
with the theoretical predictions given by Eq. 7b and multi¬ 
plied by the modification factor, K, given by Eq. 16. The 
following material properties were used: E = 7880 MPa and 
( g f = 0.16 N/mm as determined by testing, G = £718 = 
438MPa, and/, = 1.2MPa. The perpendicular-to-grain ten¬ 
sile strength is highly volume dependant, and direct deter¬ 
mination by testing is questionable because an appropriate 
size and shape of the test specimen is not obvious, f has 
therefore not been determined by testing, but used as a 
(constant) “fitting” parameter. Using f = 1.2MPa, the 
above given failure loads of the plate joint tests lead to PJ 
bhj t = 3.42 for h e = 20mm, and PJbhj t = 2.39 for h e = 
40mm. These values are also indicated in Figs. 8-11, and as 
expected, the failure loads of the plate joints are very close 
to those of the long beams for small values of a, while the 
plate joints give conservative results for large values of a. As 
shown in Fig. 6 , a constant overhang (distance from support 
to beam end) of 50mm was used in the beam tests. For small 
values of L//i e , however, it was suspected that the overhang 
had a significant influence on the failure load (especially for 
small a values). Supplementary tests were therefore con¬ 
ducted on beams with h = 200 mm and /.//z c = 2.5 for both 
h e = 20mm and h c = 40mm, using an overhang of just 
10 mm. The test results are shown in Fig. 8 a,c, which shows 
the significant influence of the overhang. 


Failure mechanism 

A majority of the long beams (L//; e = 5) failed suddenly 
without significant stable crack propagation. In a few speci¬ 
mens a crack appeared and propagated at increasing load to 
a certain critical length, after which sudden failure occurred. 
The cracks did not extend to the vicinity of the supports. 

For short beams ( Llh e = 2.5) stable crack growth was 
observed in all cases. A crack appeared and propagated all 
the way to the supports at still increasing load. After having 
reached the supports, the crack seemed to be arrested at (in 
some cases) still increasing load, after which a sudden de¬ 
crease in load was observed as the crack propagated past 
the support to the end of the beam. The failure loads given 
in Fig. 8 are the maximum loads obtained before this sud¬ 
den decrease in load. After propagation of the crack to the 
beam end, the load in some cases increased again to finally 
end in a bending failure of the beam with depth h e at a load 
level slightly higher than shown in the figures; in other cases 
bending failure occurred when the crack suddenly propa¬ 
gated past the supports to the beam ends. It is obvious that 
the overhang of 50 mm affects the test results for specimens 
with L//? e = 2.5. With 10-mm overhang, there was only little 
increase in load after the crack had reached the supports. 

It was expected that the tests would show stronger influ¬ 
ence of the supports than predicted by the model, because 
the model only includes beam actions; arch effects were 
expected for relatively small L//z e values. However, the 
tests, in agreement with the model, show no signs of influ¬ 
ence of the supports (when eliminating the effect of over¬ 
hang) down to Llh c = 2.5. 

Tests were also conducted on specimens with L//z e = 
1.25. However, these tests were abandoned because exces¬ 
sive embedment occurred without crack propagation. This 
indicates that the splitting strength, in agreement with the 
theoretical predictions, increases significantly for L//? e in the 
range between 1.25 and 2.5. 


Conclusions 

A quasi-non-linear fracture mechanics model based on 
beam on elastic foundation theory was applied to the analy¬ 
sis of splitting failure in simply supported beams with 
dowel-type fastener connections loaded perpendicular to 
grain. Tests were conducted on beams with a single dowel at 
midspan. Of particular interest was the question: how long 
beams are needed to eliminate the influence of the supports 
on the failure load? 

Tests as well as theoretical predictions show that the 
supports have no influence for span/edge distance ratios 
(2L//z c ) above 5. The theory predicts that for span/edge 
distance ratios of less than about 4, the supports have a 
strong influence on the failure load. For a span/edge dis¬ 
tance ratio (2 Llh e ) of 2.5, tests resulted in constantly in¬ 
creasing embedment without splitting failure. 

The tests indicate that the theoretical model agrees very 
well with experiments not only in terms of the influence of 
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Fig. 8. Failure load as a function of span for a h c = 40mm, h = 200mm; b h c = 40mm, h = 80mm; c h c = 20mm, h = 200mm; d h c = 20mm, 
h = 100 mm. For test results, mean values ± standard deviations are given 


the supports, but also with regard to the influence of edge 
distance, h e , and total beam depth, h. The model presented 
includes important previous models such at the Van der 
Put/Leijten model 2 and the Gustafsson/Larsen model 5 as 
special cases. 
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